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ABSTRACT 

We investigate molecular evolution in a star-forming core that is initially a 
hydrostatic starless core and collapses to form a low-mass protostar. The re- 
sults of a one-dimensional radiation-hydrodynamics calculation are adopted as a 
physical model of the core. We first derive radii at which CO and large organic 
species sublimate. CO sublimation in the central region starts shortly before the 
formation of the first hydrostatic core. When the protostar is born, the CO sub- 
limation radius extends to 100 AU, and the region inside < 10 AU is hotter than 
100 K, at which some large organic species evaporate. We calculate the temporal 
variation of physical parameters in infalling shells, in which the molecular evo- 
lution is solved using an updated gas-grain chemical model to derive the spatial 
distribution of molecules in a protostellar core. The shells pass through the warm 
region of 10 — 100 K in several x 10^ yr, and fall into the central star ~ 100 yr af- 
ter they enter the region where T > 100 K. We find that large organic species are 
formed mainly via grain-surface reactions at temperatures of 20 — 40 K and then 
desorbed into the gas-phase at their sublimation temperatures. Carbon-chain 
species can be formed by a combination of gas-phase reactions and grain-surface 
reactions following the sublimation of CH4. Our model also predicts that CO2 is 
more abundant in isolated cores, while gas-phase large organic species are more 
abundant in cores embedded in ambient clouds. 
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Subject headings: stars: formation — ISM: molecules — ISM: clouds — ISM: 
individual (IRAS 16293-2422, L1527) 



Introduction 



In a last decade, great progress has bee n made in our understan d ing of the chemi 



cal e volution of low-mass star-forming cores (jPi Francesco et al. 1 120071 : iBergin fc Tafalla 



20061 . and references therein). Among the observational advances are the detection of chem- 
ical fractionation in several prestellar cores; emission lines of the rare isotopes of CO are 
weaker at the core center than at outer radii, while emission lines of nitrogen-bearing species 
e.g., N-^H"*") show relatively good correl ation with the centrally-peaked dust continuum (e.g. 



CaseUi et al. Ill999l : iTafalla et al. Il2002[ ). 



Theoretical rnodels show that the CO depletion is caused by adsorption onto grains 
(IBergin fc Langer 1119971 : 1 Aikawa et al. Il200ll ). and that the CO depletion helps to temporar- 
ily maintain the N2H"'" abundance at the core center. Theoretical models also predict that a 
fracti on of the adsorbed CO w i ll be hydrogenated to form H2CO and CII3OH on grain sur- 
faces ( Allen fc Robinson 1977: iHasegawa et al. 1992), wh ich is confirmed by laboratory ex- 



periments ( jWatanabe fc Kouchi 



2002a 



Fuchs et al. 



20071 ). The existen ce of solid methanol 



i n low -mass star formation regions has been confirmed observationally; iPontoppidan et al. 



(I2OO3I ) detected a high abundance of CH3OH ice (15-25 % relative to water ice) towards 
three low-mass protostars among ~ 40 observed protostars. R adio observations find gaseous 



CH3OII to be abundant in the central regions of proto stars (Schoier et al 



the formation of CH3OH is inefficient in the gas phase (IGarrod et al. 



2006 



2002a|). Since 
2007h . it must 



be formed by the hydrogenation of CO on grain surfaces in the prestellar core stage, and 
then sublimated to the gas-phase as the core is heated by the pr otostar. Although CH3OH 
ice is not detected towards the maj ority of low-mass protostars (IPontoppidan et al. II2003I ) 



and the background star Elias 16 (iChiar et al. 1 119961 ). the upper limit for the CH3OH ice 
abundance is a few % relative to water ice, which is not low enough to contradict the idea 
that gaseous CH3OH around protostars is originally formed by grain-surface reactions and 
then sublimated. 

Two major questions, however, remain to be answered, the first being at what stage 
CO returns to the gas phase. The core remains nearly isothermal as long as cooling by 
radiation is more efficient than heating by contraction (compression). Eventually, though, 
the heating overwhelms the cooling, so that the core center becomes warmer. A newly-born 
protostar further heats the surrounding core. Laboratory experime nts show although some 
fraction of CO can be entrapped in water ice (ICoUings et al. II2003I ). a significant amount of 
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CO sublimates at around 20 K (jSandford &: AUamandola Ill988l ). In order to predict if an 
observable amount of CO returns to the gas-phase during the prestellar core stage or after the 
birth of a protostar, the temperature distribution in a core should be calculated by detailed 
energy transfer. Such a prediction is important in order to observe very young protostars. 
Once it is sublimated, CO again becomes a useful observational probe. In addition, CO 
sublimation significantly affects the gas-phase chemistry; for example, it destroys N2II+. 

The second question concerns how large organic molecules are formed in protostellar 
cores. In recent years, diverse organic molecules, including methanol (CH3OH), dimethyl 
ether (CH3OCH3), acetonitri le (CH^CN), and forni ic acid (HCOOH), have been detected 
towards low- mass protostars (jCeccarelli et al. 2007 . and references therein). They are still 
referred to as "hot-core species", because it was once thought that they are only characteris- 
tic of hot (T ~ 200 K) cores in high-mass star forming regions. A large number of modeling 
studies on hot-core chemistry show that sublimed formaldehyde (H2CO) and CH3OH are 
transformed to ot her organic species by g as-phase reactions within a typical timescale of 
10^ — 10^ yr (e.g. iMillar fc Hatchell Ill998l ). In low-mass cores, however, the time scale for 
the cloud material to cross the hot (T ~ 100 K) region should be sr naller than 10"^ yr, consid - 
ering the infall velocity and temperature distribution in the core (IBottinelli et al. Il2004al ). 
Furthermore, theoretical calculations and laboratory experiments recently showed that gas- 
phase reactions are much less efficient in producing some hot-core species, suc h as methyl 



formate (HCOOCH3) and dimethyl ether, than as sumed in previous models (IHorn et al. 
20041 : ICeppert et aLlbood : ICarrod fc Herbst hood ). 



Several model calcul ations have b een p erformed on the chemistry that occurs in low- 
mass protostellar cores. iDoty et al. I (120041 ) solved a detailed gas-phase reaction network 
assuming a core model for IRAS 16293-2422, and succeeded in reproducing many of the 
observed lines within 50 %. The physical structure of the core; i.e., its density and tempera- 
ture distribution, was fixed with time. They assumed gas-phase i n itial m olecular abundances 
that pertain to the high- mass hot-core AFGL 2591. iLee et al. I (120041 ). on the other hand, 
constructed a core model that evolves from a cold hydrostatic sphere to a protostellar co re by 



combi ning a sequence of Bonnor-Ebert spheres with the inside-out collapse model by IShu 
( 119771 ). They solved a chemical reaction network that includes gas-phase reactions along 



with adsorption and desorption of gas-phase/ice- mantle species. The resulting molecular 
distributions and line profiles are significantly different f r om th ose of the static core models 
( iLee et al. I l2005l ). The chemical network of iLee et al. I (120041 ). however, does not include 
large organic species. 

In the present paper, we re-investigate molecular evolution in star-forming cores with 
the partial goal of answering the two questions posed above. We adopt a core model by 
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Masunaga &: Inutsuka I (120001 ): it accurately calculates the radial distribution of tempera- 
ture, which determines when and where the ice components are sublimated. The model 
also enables us to follow molecular evolution smoothly from a prestellar core to a proto- 
st ellar core. The chemistry includes both gas-phase and grain-surface reactions according 
to iGarrod fc Herbst I (120061 ): the surface reactions in particular are important for producing 



organic molecules in a warming environment. Here, we report a solution of the reaction 
network following the temporal variation of density and temperature in infalling shells to 
derive a spatial distribution of molecules, including complex organic ones, in a protostellar 
core. 



2. Model 

2.1. Physical evolution of a star-forming core 

Figure [1] schematically shows the evolution of a star-forming core. As a model for 



this p rocess, we adopted and partially reran the model calculated by lMasunaga fc Inutsuka 



(120001 ). These authors solved the non-gray radiation hydrodynamics to follow the core evolu- 
tion from a dense starless (prestellar) core to a protostellar core assuming spherical symmetry. 
Conservation equations of mass, momentum, and energy were coupled with the frequency- 
dependent radiation tr a nsfer, which was solved by the variable Eddington factor method. 



Masunaga &: Inutsuka I (120001 ) calculated two models with different initial conditions: a ho- 
mogeneous core (one of uniform density) and a hydrostatic core (a Bonnor-Ebert sphere). 
We chose the latter one in the present work. 

Initially, the central density of the hydrostatic prestellar core is 1.415 x 10^^^ g cm~'^, 
which corresponds to a number density of hydrogen nuclei nn ~ 6 x 10'^ cm"^. The outer 
boundary is fixed at r = 4 x 10^ AU, so that the total mass is 3.852M0, which exceeds the 
critical mass for gravitational instability. The temperature in the core is initially ~ 7 K at 
the center and ~ 8 K at the outer edge; the cooling by dust thermal emission is balanced 
with the heating by cosmic rays, cosmic background radiation, and ambient stellar radiation. 



In the original model by iMasunaga &: Inutsuka I (120001 ) . the core starts contraction im- 



mediately. In the present work, however, we assume that the core keeps its hydrostatic 
structure for 1 x 10^ yr, implicitly assuming that turbulence supports it. This period sets 
the initial molecular abundances for the collapse stage. After 1 x 10^ yr, the core starts to 
contract. The contraction is almost isothermal as long as the cooling rate overwhelms the 
compressional heating, but eventually the latter dominates and raises the temperature in the 
central region. Increasing gas pressure decelerates the contraction and eventually makes the 
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at the center. When the core center becomes 



7 g 3 ^-^^ 2000 K, the hydrostatic core becomes unstable due 



first hydrostatic core, known as "the first core 
as dense as 10 

to H2 dissociation and starts to collapse again (the second collapse). The central density 
increases rapidly, and within a short period of time the dissociation degree approaches unity 
at the center. Then the second collapse ceases and the second hydrostatic core; i.e. the 
protostar, is formed. The protostar is surrounded by the infalling envelope, which we call 
the protostellar core. After the onset of contraction, the initial prestellar core evolves to 
the protostellar core in 2.5 x 10^ yr. After the birth of the protostar, the model further 
follows the evolution for 9.3 x 10'^ yr, during which the protostar grows by mass accretion 
from the envelope. At each evolutionary stage, the model gives the total luminosity of the 
core and the radial distribution of density, temperature, and ii ifall velocity at r > 10~^ AU. 
More detailed explanat i ons o f core evolution can be found in iMasunaga et al. I (Il998l ) and 
Masunaga &: Inutsuka I (120001 ). 



In the present paper, we report a re-analysis of their results in order to find the sublima- 
tion radii of CO and large organic species at e ach ev olutionary stage. We obtained electronic 
data and the code of IMasunaga fc Inutsuka I (I2OOOI ) from the author (H. Masunaga), and re- 
ran the calculation of the prestellar core phase to increase the number of time steps at which 
the core structure is recorded. As for the protostellar phase, we used the existing data. We 
are interested in the regio n outside 1 AU, which is of imp ortance for radio observations but is 
not discussed in detail by IMasunaga &: Inutsuka I (120001 ). Based on this physical core model. 



we followed the temporal variation of density, temperature, and visual extinction in in falling 
shells, in which the chemical reaction network was solved (lAikawa et al. Il200ll . l2005l ). 



2.2. Chemical reaction network 



The chemical ev olution was computed w ith the gas-grain model developed in the OSU 
astrochemical group (jGarrod fc Herbst II2006I ). The current version of the model calculates 
the abundance of 655 species (458 gas-phase and 197 grain-surface atoms and molecules) 
through a total of 6309 reactions. The model follow s the gas-phase as well a s the grain-surface 



chern istry using the rate equation method (e.g. iHasegawa et al. I Il992l : iRuffle &: Herbst 



2OOOI ) . The gas-phase network is based on the osu.2005 database available online (http:/ /www.physics.ohio- 
but without the element fluorine, which has been recently added to osu.2005. The model in- 
cludes chemical reactions that can be important at T < 300 K. At the final stage of our core 
model, the temperature is somewhat higher than this limit at r < 10 AU. In the following, 
we mainly discuss the molecular evolution that occurs at T < 300 K. 



The assumed elemental abundances, known as "low-metal" values because of their strong 
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depletions for elements heavier than oxygen, are listed in Table [H Initially the species are 
assumed to be in the form of atoms or atomic ions except for hydrogen, which is entirely in 
its molecular form. The cosmic-ray ionization rate ( is set to 1.3 x 10"^^ s~^. 

The neutral species in the gas-phase stick to the grains upon collision with a probability 
of 0.5. At high temperatures, this probability is much too high. However, since evapo- 
ration is so rapid at such high temperatures, it is unimportant what assumption is made. 
The grain-surface species can evaporate through thermal evaporation and two non-thermal 
desorption processes. The biriding e nergy of each molecular species on the grain surface is 



taken from iGarrod &: Herbst 



(120061 ). The first non-thermal mechanism is the action of a 
cosmic ray encountering a grai n, which rapidly increases the grain temperature and allows 
the evaporation of the specie s (IHasegawa &: Herbst 1119931) . The second process, as recently 
developed quantitatively by Garrod et al. (2006^2005^ is based on the assumption that 
the energy released by exothermic association rea ctions on grain surfaces can allow the par- 
tial evaporation of the products (jWilliams 1 119681 ). In this process, the fraction of product 
evaporation is given by 



/ 



a(l — ED/EreacY 



1 + a(l - Eo/EreacY-^ ' 

where E^ is the binding energy of the product, -Ereac is the energy of formation in the 
reaction, and s is the number of vibrational modes in the molecule/surface-bond system. 
The parameter a is the ratio between the surface-molecule bond frequency and the frequency 
at which the energy is lo st to the surface. We assumed an a value of 0.01 for all species based 
on lGarrod et al. I (120071 ). With this value of a, less than 1% of \ yater is evaporated du r ing its 



formation on grains. Using molecular dynamics calculations, iKroes fc Andersson I (120061 ) 



showed that it is unlikely that this fraction is larger than 1%. If it were, however, only the 
gas-phase abundances of some large molecules such as CH3OH would be increased in the 
outer parts of the protostellar envelope. 

Grain-surface species are dissociated by penetrating UV radiation and cosmic-ray in- 
duced UV radiation. We note that our model indirectly includes photodesorption via the 
second non-thermal desorption process discussed above. For example, the photodissociation 
of water ice produces OH radicals, a fraction of which recombine with hydrogen atoms to 
re-generate water. Since we assume the partial desorption of the surface-reaction products, 
UV irradiation of water ice results in desorption of water. 

Species can diffuse on the grain surfaces by thermal hopping and react with each other 
when they meet. No quantum t unneling was assuni ed even for hydrogen atoms based on 
recent experimental results (see iGarrod et al. 1 120071. for discussion on this point) . Other 
details concerning the model can be found in IGarrod fc Herbst I (120061 ) and 
2007h . 



Garrod et al. 
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During the first 1 x 10® yrs of the integration, we computed the chemical abundances 
under static dense cloud conditions in order to obtain the initial abundances for the collapse. 
The chemical evolution in infalling shells was then computed as the shells experience tempo- 
ra l variations of density, t e mper ature and visual extinction A^. Although the original model 



of iMasunaga fc Inutsuka I (120001 ) gives the temperatures of gas and dust of various composi- 
tions separately, we assumed that the dust temperature is equal to the gas temperature for 
simplicity. 

We calculated the column density of hydrogen nuclei (A^h) from the core outer edge to 
each shell and converted it to by the formula A^ = Ah/(1-59 x lO^^cm"^) mag. Initially 
the visual extinction from the outer edge (r = 4 x 10^ AU) to the core center is about 5.5 
mag, and the shells we follow are located at ~ 1 x 10^ AU, where A^ ^ 1 mag. Assuming 
that our model core is embedded in ambient clouds, we added 3 mag to the visual extinction 
obtained above, and ignored photodissociation of CO and H2. In §4, we also report molecular 
evolution in an isolated core, in which the extra 3 mag is not added to A^ 



3. Results 

3.1. Physical evolution of the core and sublimation radius 

Figure [2] shows the radial distributions of density, temperature and infall velocity in the 
model core at assorted evolutionary stages. We define the moment of the protostar formation 
as tcore = 0. The left panels refer to the prestellar phase, which occurs before the birth of 
the protostar, and the right panels to the protostellar phase. The distributions are shown at 
different times in each panel. 

We concentrate on the region outside 1 AU, which is of impor tance for radio observations 



but is not discussed in detail by IMasunaga fc Inutsuka I (I2OOOI ). At tcore < — 1 x 10'^ yr, 
more than 1 x 10'^ yr before the formation of a protostar, the core is almost isothermal 
and the core contraction is similar to the Larson-Penston (L-P) collapse with a constant 
temperature ( Larson IIiQGqI ). At tcore ~ — 1 x 10^ yr, the core center starts to heat-up, and 



the core deviates from the L-P core. The increasing pressure decelerates the contraction, 
and the first hydrostatic core, of radius ~ 1 AU, is formed at tcore ~ —5.6 x 10^ yr (see the 
sharp velocity gradient at ~ 1 AU in Figure [2]). When the central density and the central 
temperature reach nn ~ several 10^^ cm~^ and ~ 2000 K, dissociation of H2 leads to the 
second collapse, which ceases in a short time scale (~ 1 yr). In the protostellar stage, the 
density in the core envelope decreases with time because of the accretion to the central star, 
while the temperature and infall velocity increase. 
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In Figure [3l the temperature at the core center through 200 K is plotted as a function of 
the central density n^"*''^^. The dashed line indicates a temperature of 20 K, which is roughly 
the sublimation temperature of CO, or the temperature at which the rate coefficients of 
accretion and thermal desorption for pure CO ice becomes equal. In the prestellar phase, 
the core temperature is determined by the balance between cooling by thermal radiation 
and heating by cosmic rays, cosmic background radiation, ambient stellar radiation, and 
compressional effects. The central temperature for a time decreases to ~ 5 K, because the 
increasing column density of the cloud prevents the penetration of optical photons. The 
core temperature starts to rise at > 10^ cm^^, when compressional heating dominates 

radiative cooling. The temperature increase is, however, rather moderate until the central 
density reaches ~ 10^^ cm~^. Shortly thereafter, CO sublimation starts when the central 
density reaches a few x 10^^ cm~^. This occurs only several x 10^ yr before the first-core 
formation and ~ 10^ yr before the second-core formation. 

Figure H] shows the temporal variation of two radii, inside of which the temperature is 
higher than 20 K (r20K) and 100 K (rioox)- The former is roughly the sublimation temper- 
atures of CO. Desorption energies and thus sublimation temperatures vary among species, 
but some large organic species such as CII3OII and HCOOH, which are of current interest, 
have sublimation temperatures of ~ 100 K. The sublimation radii are plotted as a function 
of total luminosity of the core instead of the core age, since the former is an observable value. 
The labels with arrows depict the core age tcore- The plot starts about 770 yr before the 
second-core formation, when the first core is about to be formed. As the first core grows 
in density, the core luminosity and envelope temperatures increase. At the moment of the 
second-core formation (tcore = 0), r20K and riooK are ~ 100 AU and ~ 10 AU, respectively. 
At the final stage of the model, tcore = 9.3 x 10^ yr, the total luminosity reaches ~ 23L0, 
and the sublimation radii are r20K ~ 3.9 x 10^ AU and riooK ~ 1-2 x 10^ AU. 



3.2. Physical conditions in infalling material 



Up to now, we have described the core model with physical parameters that are given 
as functions of radius at each evolutionary stage, the so-called Eulerian approach. In order 
to calculate molecular evolution, however, it is more convenient to follow the Lagrangian 
approach, in which we consider the temporal variation of physical parameters in assorted 
collapsing shells, because this tem poral variation represents the conditions in which the 
chemical reaction network is solved (lAikawa et al. II2OOII . l2005l ). 



Figure [5] shows the temporal variation of density and temperature in shells that reach 
r = 2.5, 15, 125 and 500 AU at the final stage of our model (tcore = 9.3 x 10^ yr = tgnai)- 
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The infalling shells begin at r ~ 1 x 10^ AU in the prestellar core, where they stay for 
1 X 10^ yr, a temporal period omitted in Figure [51 After the core starts contraction, each 
shell migrates inwards. In the prestellar stage, the variation of density and temperature is 
not very significant; these shells are at r > 6 x 10^ AU, where the infall velocity is small and 
the radial gradient of density and temperature is moderate (Figure [2j). In the protostellar 
stage, which occupies a large fraction of Figure [5l they migrate to inner radii where the infall 
velocity, temperature, and density steeply increase inwards. Hence, the temporal variation 
of density and temperature accelerate. In order to highlight the rapid variation near and at 
the final stage, the horizontal axis in Figure [5] is set to be the logarithm of tfinai — ^core- 

A key parameter for the chemistry of hot-core species is the time scale of the hot-core 
phase and/or, its predecessor, the warm-up phase. In a classical hot-core chemical model, in 
which the gas-phase reactions of sublimates are followed with a fixed temperature of ~ 200 
K, it takes 10^ — 10^ yr to form new large organic species. In other words, the hot-core 
phase of 10^ — 10^ yr is needed to account for the existence of large organic species (e.g. 



Millar fc Hatchell Ill998l ). In more recent hot-core models, the temperature first increases 
with time rather than remaining fixed; the warm-up phase is the period in which the tem- 
perature rises from 10 K to ~ 100 K. Various reactions occur during the warm-up phase, 
because many volatile species ret urn to the gas phase, and because the sublimation tem- 
perature varies with species (e.g. IViti &: Williams Ill999l ). Grain-surface reactions between 
molecular radicals conta ining heavy elements are also enhanced in this temperature range 



(iGarrod fc Herbst II2006I ). Considering the formation time scale of stars, IViti et al. I (120041 ) 
argued that the warm-up phase is longer in stars with lower masses, and estimated the 
warm-up phase to be 10^ yr for solar-type stars. Our model, in contrast, gives a warm-up 
phase of only several 10^ yr; since the gas and dust are falling inwards towards the central 
star, the timescale of the warm-up phase is determined not by the time scale of the star 
formation, but by the size of the warm region divided by the infall velocity. The timescale 
of the hot-core phase, in which the temperature of the infalling shells is ~ 100 — 200 K, 
is even shorter; the shell reaching r = 2.5 AU at tfinai spends only 100 yr in the region of 
T > 100 K before it falls onto the central star. There would thus not be sufficient time 
for the gas-phase reactions of sublimated species to form complex species. We are dealing 
here with the formation of a low-mass protostar, so the hot-core stage is more aptly referred 
to as a "hot coririo", and short time scales hav e previously been discussed for such regions 
(IWakelam et al. I |2004J : iBottinelli et al. Il2004al ). It should also be noted that in our con- 
tracting core model, the central region with high temperatures (T > several tens of K) is 
continuously fed by infalling materials, which are then lost to the central star. 
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3.3. Molecular evolution in infalling shells 

Figure [6] shows the temporal variation of molecular abundances in the infalling shell 
that reaches r = 2.5 AU at tgnai- The horizontal axis is again the logarithm of tgnai — ^core- 
Black and gray lines represent species in the ice mantle and in the gas phase, respectively. 
The chemical evolution at tgnai — ^corc < 4 yr should be taken with caution because the 
temperature in the shell exceeds 300 K, while our chemical network is originally constructed 
for temperatures of < 300 K. The species shown in Figure [6] are mostly (except CO and 
H) formed by grain-surface reactions. They are desorbed into the gas phase when the shell 
temperature reaches the sublimation temperature of each species; at this time the fractional 
abundance of the surface species is simply transferred to a similar fractional abundance in 
the gas. 

Water, the most abundant molecule with a heavy atom, is formed by the hydrogenation 
of oxygen atoms on a grain surface, and is already abundant when the contraction starts. 
Methanol, on the other hand, takes longer to reach a high abundance: it is formed by 
the hydrogenation of solid CO, which occurs efficiently until the shell temperature reaches 
> 15 K and the H atom abundance in the ice mantle drops due to evaporation. The 
fractional abundance of methanol exceeds 10^^ at times within 10^ yr of the final time. 
Then the abundance reaches a constant value; CH3OH is constantly destroyed by the cosmic- 
ray induced photodissociation and re-formed by a low rate of hydrogenation. Solid carbon 
dioxide is formed by O + HCO ^ CO2 + H and CO + OH ^ CO2 + H. It reaches a 
fractional abundance of 10~^ about 10^ yr before the final time. 

The large organic species are also formed mainly by grain-surface reactions. Dimethyl 
ether (CH30CH3)and methyl formate (HCOOCH3) are produced by the surface reactions 
of the radical H2COH with CH3 and HCO, respectively, where H2COH is formed by hy- 
drogenation of H2CO. Formaldehyde, at this stage, is constantly formed by hydrogenation 
and destroyed by a reaction with OH ice and cosmic-ray induced photodissociation. Both 
CH3OCH3 and HCOOCH3 reach terminal fractional abundances of ~ 10~^. Formic acid 
(HCOOH), on the other hand, is formed by the gas-phase reaction of OH with sublimated 
H2CO, after which it is adsorbed onto grains, and remains in the ice mantle until the temper- 
ature gets high enough for its sublimation. Its asymptotic fractional abundance is slightly 
greater than 10^®. Comparing Figure [H] with Figure we can see that these large organic 
species are mostly produced at temperatures of 20 — 40 K, rather than at T > 100 K where 



CH3OH is sublimated into the gas phase (see also iGarrod fc Herbst 1120061 ). In this tem- 
perature range, heavy-element species such as CO, HCO and CH3 can diffuse on the grain 
surface and form complex molecules efficiently. 
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3.4. Distribution of molecules in a protostellar core 

In order to derive the spatial distribution of molecules in the protostellar core, we have 
calculated molecular evolution in 13 shells, which reach r = 2.5 — 8000 AU at tfinai- Figure [7] 
shows the radial distribution of (a) physical parameters and (b-d) molecular abundances at 
^finai- Abundances inside 10 AU are not shown for the following two reasons: firstly they are 
very similar to the abundances at 10 AU, except that the H2CO abundance drops inwards, 
and secondly, the temperature in the inner radius is somewhat higher than the temperature 
range originally considered in our chemical model (§2.2). Figure [7| (b-d) show that large 
organic molecules are abundant in the gas phase in the central region of the protostellar core. 
Since they are formed mostly on grain- surfaces, rather than from gas-phase reactions among 
sublimates, their gas-phase abundances sharply increase at radii corresponding to their own 
sublimation temperatures and remain high at smaller radii. For example, CH3CHO extends 
to ~ 500 AU, while CH3OH extends to ~ 100 AU. One exception, again, is HCOOH, which 
extends to ~ 500 AU because of its formation in the gas phase. 

Large organic species are abundant in ice mantles at radii of 100 — 1000 AU. For radii 
r > 10^ AU, the abundances of these species tend to decrease sharply with increasing radius 
except for methanol. Water, another hydrogen-rich, or saturated, species, remains nearly 
constant in abundance. In the outer re gion (r > 4000 AU), a large fraction of carbon on the 



grain surface is in the form of CH4 (cf iGarrod et al. 1120071 ). 



4. Discussion 



4.1. Physical structure of the core and sublimation radius 



We have re-analyzed and adopted the results of iMasunaga fc Inutsuka I (I2OOOI ) to derive 
the sublimation temperatures of CO and large organic species, and to investigate molecular 
evolution in a star-forming core. The chosen conditions in the envelope can be different from 
those of other models because the temperature distribution in the envelope depends not only 
on the evolutionary state and mass of the central object (either the first or second core), but 
also on t he mass distribu tion in the envelope, which should depend on the initial conditions. 
Recently, lOmukai I (120071 ) investigated the temper ature distribut ion in the first-co re envelope 
assuming mass distributions from the L-P model (ILarson Ill969l ) and Shu model (IShu 1119771 ) 
similarity solutions. The former has a more massive envelope than the latter. When the 
first core has a mass of 0.05Mq, for example, the hydrogen number density at r = 10 AU is 
several 10^° cm~^ and ~ 10^ cm~^ with the L-P model and the Shu model, respectively. The 
L-P model gives a core luminosity of 1O~^L0, and sublimation radii of r20K ~ 100 AU and 
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'"looK ~ 10 AU, while the Shu model yields a core luminosity of ~ 1O"^L0 and an r20K of 
~ 20 AU. The latter model does not exceed 100 K at any radius. Our core model is warmer 
than the Shu model but colder than the L-P model. 

In the second core stage, on the other hand, the envelope structure can deviate from 
spherical symmetry and be accompanied by a circumstell ar disk. Considering a typica l an- 



gular velocity for molecular cloud cores (~ 10"^^ s~^; e.g.. iGoldsmith fc Arquilla Ill985l ). the 
centrifugal radius (i.e. the initial disk radius) is ~ 100 AU. Our results at r > several hun- 
dred AU are thus relatively robust, while the core structure would be significantly different 
from our model at smaller radii. For example, the large organic species could be sublimated 
by the accretion shock onto the forming disk rather than in the envelope, since the centrifu- 
gal radius ~ 100 AU coincides with the sublimation radius of large organic species in our 
model. 



4.2. Simple molecules 



Lee et al. I (120041 ) investigated the evolution of relatively simple molecules, such as HCN 
and N2H"'", in a star-forming core by combining a se quenc e of Bo nnor-Ebert spheres for the 
prestellar stage with the inside-out collapse model of IShu I (119771 ) for the core after the first- 
core formation. These simple molecules are often observed in star-forming cores and are of 
importance as observational probes, while we mainly discussed large organic species in §3. 
Figure [8] shows the radial distribution of simple molecules in our model at tsnai- 



Comparison with iLee et al. I (120041 ) is not easy because there are many differences in 
the physical core mo dels and chemica l reaction networks. First, the inside-out core model, 
which is adopted in iLee et al. I (120041 ). yields lower temperatures than our core model, as 



discussed above. Secondly, we adopt different binding energies for molecules to the grain 
surface. While the difference is on the order of only a few hundred K for many species, the 



binding energy of NH3 is significantly higher in ou r model (^5534 K) than inlLee et al. I (120041 ) 



(1082 K); the latter value seems to originate from iHasegawa fc Herbst I (119931 ) . in which the 
hydrogen bonding of NH3 is not taken into account. Thirdly, grain-surface reactions, which 
are the main forma tion processes of s aturated species such as NH3 and H2CO in our model. 



are no t included in iLee et al. I (120041 ). Hence, here we compare our results with iLee et al. 
(120041 ) qualitatively rather than quantitatively. 



A main conclusion of Lee et al 



(120041 ) is that the molecular abundances vary signif- 
icantly near and inside the CO sublimation radius. For instance, N2H"'" is destroyed by 
CO, and thus declines steeply inwards across the CO sublimation radius, an effect which 
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happens in our model as well. iLee et al. I (120041 ) also found that some molecules, such as 
H2CO, have their peak abundance at the sublimation radius; the gas-phase abundance first 
increases inward via sublimation, but then decreases due to gas-phase reactions at smaller 
radii. Similar phenomena can be seen i n our model , but th e spatial variation of the gas-phase 



abundances is more moderate t han in iLee et al. I (120041 ). for which we can think of a few 



reasons. First, iLee et al. I (120041 ) calculated the molecular evolution in 512 shells, while we 
calculated the chemistry in only 13 shells. A larger number of shells are needed to resolve 
abundance fiuctuations on a smaller radial scale. Secondly, higher infall speeds, in general, 
make the abundance distributions more uniform. Since the infall speed is higher at inner 
radii in the infalling envelope, it would be natural that th e molecular abundances remain 
relatively uniform at r < 100 AU, which is not calculated in iLee et al. I (l2[)04f). Third l y, our 



chemical network includes a much larger number of species and reactions. ILee et al. I (120041 ) 
used a reduced reaction network with ~ 80 species and ~ 800 reactions to save computa- 
tional time, while our model includes 655 species and 6309 reactions. In a small reaction 
network, a sudden increase of a species (e.g., as caused by sublimation) can easily change the 
abundances of other species though chemical reactions. In a large reaction network, on the 
other hand, a larger number of reactions contribute to the formation and destruction of each 
species, and thus the sudden abundance change of one species does not necessar arily prop- 



agate to other species. Another noticeable difference in our results from those of iLee et al. 



(I2OO4I ) is that HCO"'" decreases more steeply inwards at r < 1000 AU in our model; H3CO"'" 
and CsH^ are the dominant positive ions rather than HCO^ at the inner radii. 

The oxygen-bearing species atomic oxygen (O), molecular oxygen (O2), and gaseous H2O 
are also of special interest because O is very reactive and because the two molecules have 
been intensively observed by the Suhmillimeter Wave Astronomy Satellite (SWAS) and Odin 
Satellite in recent years. Our model predicts that O reaches an abundance of 4 x 10~^ relative 
to hydrogen nuclei at r = 8000 AU, while it steeply decreases inwards from ~ 1000 AU to 
100 AU. E ven at r = 8000 AU , however, H2O ice is the most abundant 0-bearing species 
(Figure!?]). iBergin et al. I (l2000l ) summarized the 5* W/IS* observations of cold molecular clouds 
by stating that the fractional abundance of O2 lies under 10~^ and that of gaseous water lies 
in the range 10~^ to a few x 10^^. Molecular oxygen has recently b een detected toward s 
p Oph by Odin with an abundance of 5 x 10~^ relative to hydrogen (ILarson et al. II2007I ). 



These abundances are consistent with our predictions for the outermost radius r = 8000 



AU, where the density nn 
molecular clouds. 



10^ cm ^ and temperature T < 20 K are similar to those in 
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4.3. Dependence on visual extinction at the core edge 



So far we have assumed that the model core is embedded in ambient clouds of = 3 
mag. In reality, some cores (e.g. Bok globules) are isolated, while others are embedded 
in clouds. Isolated cores are directly irradiated by interstellar UV radiation, which causes 
photo-reactions (photodissociation and ionization) in the gas-phase and ice mantles (e.g. 
Lee et al. I EoOgI : iRuffle fc Herbst I l200ll ) . In order to evaluate the effect of ambient U V 
radiation, we recalculated molecular abundances in a core that is directly irradiated by 
interstellar UV radiation; i.e. = mag at the outer edge of the core (r = 4 x 1 "^ AU) . 



The photodissociation rates of H2 and CO were calculated by following iLee et al. I (119961 ) . 
which gives the shielding factors as a function of A^ and column densities of CO and H2 in 
the outer radii. The outermost shell for which we calculate molecular evolution is initially 
located at r ~ 1.4 x 10^ AU and declines to r = 8000 AU at tfinai- Column densities of 
CO and II2 outside of this shell were estimated by assuming n{CO)/nu = 5 x 10~^ and 
n(H2)/nH = 0.5. 

Figure [9] shows the resultant distribution of molecular abundances in a protostellar core 
at tfinai- Compared with the embedded core model (Figure [7]), the fractional abundances of 
CH3OII and II2CO are lower by more than one order of magnitude, while the CO2 abundance 
is higher. When the shells are still at r > several thousand AU and have relatively low A^ 
(< 4 mag), the photodissociation of H2O ice efficiently produces OH, which reacts with 
CO to produce CO2 ice. Methanol in the ice mantle is dissociated to H2CO, which is 
further dissociated to CO. Species such as CH3CHO, HCOOCH3 and CH3OCH3 are also 
less abundant in the isolated model, because their formation path includes H2CO in the ice 
mantle. Formic acid in the gas phase extends only up to ~ 100 AU, while it extends to 
several hundred AU in the embedded model. In the isolated model, it is formed mainly by 
OH -|- HCO on the grain surface. Our results may indicate that molecular abundances in 
hot cores and corinos depend on whether the core is isolated or embedded in clouds. 

The importance of photo-reactions on ice mantle ab undances has also been investi gated 
in a number of laborat ory experiments. For example, IWatanabe fc Kouchi I (j2002bl ) and 
Watanabe et al. I (120071 ) found that carbon dioxide is efficiently produced by UV irradiation 
on a binary ice na i xture of H2O and CO, a result that is consistent with our model. However, 
Watanabe et al. I (120071 ) found that the UV irradiation also produces CH3OH with a relative 
abundance of ?T,(CH30H)/r2(CO) ~ 10 % in the ice mixture. Comparison of our model with 
their experiment is not straightforward because of differences in temperature and included 
reactions, but we may have underestimated the CH3OH ice abundance in the irradiated 
core model. The discrepancy can arise from the H atom desorption rate in our model. We 
calculated that UV radiation penetrates into the ice mantle and dissociates H2O to produce 
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H atoms embedded in ice. Although we do not discriminate between H atoms on the ice 
surface and those embedded in the ice mantle, the embedded H atoms in reality would have 
a lower desorption rate and a better chance of reacting with neighboring CO, a reaction 
that initiates the formation of methanol in the mantles. Discrimination between surface and 
embedded hydrogen atoms should be included in future work. 



4.4. Comparison with observation 

Our model results can be compared with observational results of low-mass protostars. 
Comparison of the physical s tructure of the core has already been discussed in detail by 



Masunaga fc Inutsuka I (120001 ). Here we concentrate on molecular abundances. First, we 



compare molecules in ice ma ntles. The observation of ices towards the low-mass proto 



star Ehas 29 is summarized in lEhrenfreund fc Shutte I (120001 ): the abundances of CO, CO2 



CH3OH and CH4 r elative to water ice a re 5.6 %, 22 %, < 4 %, and < 1.6 %, respectively. 



On the other hand, iPontoppidan et al. I (120031 ) detected a high abundance (15-25 % relative 



to water ice) of CH3OH ice towards 3 low-mass protostars among 40 observed protostars. 
Although the observation of ice features is difficult, it is probable that the composition of 
ice mantles depends on their environment and the history of the observed regions. From a 
theoretical point of view, the abundances of molecules on grains and their fractional abun- 
dances compared with H2O ice depend on time and radius from the protostar. In addition, 
an embedded core and an isolated core lead to significantly different calculated abundances 
for sohd CH3OH, H2CO and CO2. 

Table [2] lists ratios of ice species with respect to water ice at t^nai for local abundances 
at radii of 1000 AU and 8000 AU and for column densities towards the core center. Since 
most of the core material exists at r < 8000 AU along the line of sight, and since the ices are 
present at r > 10 AU, the column density is calculated by integrating the number density of 
ice species from 10 AU to 8000 AU. It is interesting to note that regardless of the distance 
from the protostar, the isolated core model gives smaller surface abundances of CH3OH and 
H2CO and a higher surface abundance of CO2 than the embedded model. In both models, 
the surface abundance of CO at 1000 AU is much smaller than observed, but it increases 
towards larger radii up to 10 and 20% in the embedded and isolated models respectively. 
Comparison with the observations is best done using our column density ratios, which are 
in reasonable agreement with Elias 29 for both models. Considering the variation among 
cores, our models show reasonable agreement with observation. The disagreement with CO 
doubtless results from our assumption concerning desorption rates. In the present work, 
we use a desorption energy for each species mainly referring to laboratory experiments on 
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pure ice sublimation or theoretical estimates that sum up the van der Waals forces between 
adsorbed atoms and grain surface. But in reality, interstellar ice is a mixture, and hence 
volatile species can be entrapped by less volatile species; recent temperature-programmed 
desorption results sho w that in a mixture r ich in water ice, much CO desorbs at considerably 
higher temperatures (jCoUings et al. 1 120031 1 . 



Table[3]lists the gas-phase abundances of large molecules towards the low-mass protostar 
IRAS 16293-2422 and in the central region (r = 30.6 AU) of our core mo del. The physical 



and ch er nical structures of IRAS 16293-2422 have been intensively studied (jCeccarelli et al. 
2000al lbl: ISchoier et ani2002bl : ICauzax et ani2003l : Ichandler etahlboosh . The physical pa- 
rameters derived by the model at tgnai (~ 23L0 and riooK ~ 1-2 x 10^ AU) ar e very close to 



the o nes of IRAS16293-2422; the observed luminosity of the source is 27 Lq (jWalker et al. 



19861 ) and the physical structure has been constrained by mu lti-line analysis of H2O and 
H2CO observations through a detailed radiative transfer code (jCeccarelli et al. Il2000bl ). It 
should be noted that the emission lines of large organic molecules observed in this source 
are not spatially resolved, except for a few lines investigated by inter ferometric observations 
jBottinelh et al. Il2004bl : Ixuan et al.lbooi iRemijian fc HoUis II2OO6I ). The density and tem- 
perature of the gas should vary both within the beam and along the line of sight. Hence the 
derived molecular abundances vary significantly depending on the assumptions concerning 
the physical structure of the core and the emitting region, which explains the difference in 
abundances determined by different investigators (see Table |3]). 

Although it is not obvious which core model, embedded or isolated, should be compared 
with IRAS16293-2422, we would prefer the embedded core model. While IRAS 16293-2422 is 
in the Ophiuchus molecular cloud, which harbors several UV sources in the for m of OB stars , 



the ^^CO and C^^O obs ervations indicate that the core is embedded in dense gas ( iLoren Ill989 



Tachihara et al. II2OOOI ). The high molecular D/ H ratios observed towards IRAS 16293-2422 
and its neighboring starless core IRAS16293E (jCeccarelli et al. 1 120071 : IVastel et al. 1 12004| ) 
indicate that these cores have been very cold and thus well-shielded from the ambient stellar 
radiation. Considering the uncertainties in observationally-estimated molecular abundances, 
HCOOCH3, HCOOH, and CH3CN in our embedded core model show reasonable agreement 
with the observations. The embedded model, however, underestimates CII3OCH3 and over- 
estimates H2CO and CH3OH. The isolated core model, on the other hand, reproduces ob- 
served abundances of CII3OH, HCOOH and CH3CN, but strongly underestimates H2CO, 
HCOOCH3, and CH3OCH3. 

We can think of several possible solutions to improve the agreement with observation. 
First, it is noteworthy that the gaseous CH3OH abundance estimated in IRAS 16293-2422 is 
much lower than the abundance of CH3OH ice (n(CH30H ice)/nH ~ 10~^, assuming n(H20 
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ice)/ riYL ~ 10~^) detected by iPontoppidan et al. I (120031 ) towards some low-mass protostars. 
We may have missed or underestimated the reactions which transform CH3OH to other large 
organic species. Secondly, the abundances of large organic species in the central region can 
vary with time. In the present work, we have concentrated on the distribution of molecules 
only at tfinai- Shells that reach the central region at different times should experience different 
temporal variations of physical conditions. Some shells may experience longer periods at 
T ~ 20 — 40 K, where large organic species start to be efficiently formed. This possibility 
will be pursued in a future work. Thirdly, core models with rotation could produce higher 
abundances of large organic species; because of the centrifugal force, core material migrates 
more slowly and stays longer in the temperature range preferable for the formation of large 
organic molecules (e.g. in a forming disk). 

Because of the beam size of the current radio observations, little is known concern- 
ing the spatial distribution of the large org anic species; they are mostly confined within a 



few arc seconds from the core center (e.g. iKuan et al. II2004I ). But there are exceptions 



Remijian fc HoUis I (120061 ) found that HCOOH and HCOOCH3 show extended emission of 
~ 5 arcsec. Our embedded core model reproduces the extended emission of HCOOH, while 
HCOOCH3 is confined to r < 100 AU. 



4.5. Carbon chains in a protostellar core 



Recently, ISakai et al. I (j2007a[ ) detected strong emission lines of carbon chain species 
such as C4H and C4H2 towards the low-mass protostar LI 527, which is considered to be in a 
transient phase from class to class I. In general, carbon-chain species are associated with 
the early stages of a cold cloud core when the dominant form of carbon changes from atomic 
carbon to CO. Hence, the existence of carbon-chain species towards L1527 is a surprise. 

Figure [10] (a-b) shows the temporal variation of CH4 and carbon-chain abundances in 
the shell that reaches r = 2.5 AU at tfinai- When the core starts to contract, methane 
(CH4) and C3H4 are already abundant in ice mantles. Methane has been formed by the 
hydrogenation of carbon on grain surfaces, while C3H4 has been formed by a combination of 
gas-phase reactions (to form unsaturated carbon chains such as C3H and C3H2) and grain- 
surface reactions (to hydrogenate them). When the CH4 sublimates, some fraction reacts 
with C^ to form C2Hj^, which is a precursor for ion- molecule reactions that lead to the 
production of larger unsaturated hydrocarbons. While gas-phase reactions tend to make the 
carbon chains longer, adsorbed species experience hydrogenation and dissociation by cosmic- 
ray induced UV radiation. Figure [10] (c-d) shows the distribution of carbon-chain species 
in our embedded core model. Methane and C3H4 are abundant in the ice mantle even at 
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r > 10^ AU, while other hydrocarbons are abundant at a few 100 AU ^ r < 10^ AU. Most 
of the carbon-chain species desorb at r < a few 100 AU. 

In summary, our model indicates that carbon-chain species can be re-generated in the 
protostellar core by the combination of gas-phase reactions and grain-surface reactions. It 
should be noted, however, that large oxygen- containing organic species are not detected in 



L1527 (jSakai et al. Il2007bl ). while they are abundant in the central region of our model. In 
future work, the temporal variation of the molecular distribution in model cores should be 
investigated to see if at certain evolutionary stages carbon-chain species are abundant while 
large organic species are not. For example, in their recent study of the chemistry of cold cores. 



Garrod et al. I (120071 ) found that generation of gas-phase hydrocarbons from precursor surface 
methane occurs at very late times, after the abundance of surface methanol has essentially 
vanished. More detailed observations are also needed for a quantitative comparison with 
models. While our model predicts the column density of C4H to be 2 x 10^^ cm~^ towards 
the central star, the observation gives 2 x 10^^ cm~^. The actual column density towards 
the central star can, however, be larger, because the emission is averaged within the beam 
(a few 10 arcsec) in the current observation. Methane ice, the precursor of the carbon-chain 
species in our model, is as abundant as 23 % relative to water ice at the outer edge of our 
core model. A deep integration of ice features towards field stars is needed to constrain the 
CH4 ice abundances at the outer edge of the core and/or quiescent clouds, while the column 
density ratio of CH4 ice to H2O ice (1 %, see Table 2) in our model is consistent with the 
observation towards YSO's (§4.3), because of the relatively large sublimation radius of CH4. 



5. Summary 



We have adopted the one-dimensional radiation-hydrodynamical core model of lMasunaga &: Inutsuka 



(I2OOOI ) to investigate molecular evolution in a low-mass star-forming core. The physical 
structure (density, temperature and infall velocity) of the core is given as a function of time 
and radius (distance from the core center). The temporal variation of density, temperature 
and visual extinction is calculated for assorted infalling shells. The temporal variation of 
molecular abundances in these infalling shells and the radial distribution of molecules in the 
protostellar core are calculated by solving a chemical reaction network that couples gas-phase 
and grain-surface chemistry. 

In the prestellar phases of star formation, species such as CO, which contain heavy 
elements, are depleted onto grain surfaces due to low temperatures in the vicinity of 10 K 
and below. As the contraction proceeds, the compressional heating overwhelms the radiative 
cooling. When the central density reaches ~ 10^^ cm~^, the temperature at the core center 
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starts to rise sharply and soon reaches the subhmation temperature of CO (~ 20 K). The 
first hydrostatic core is formed shortly (in several 10^ yr) thereafter. When the protostar is 
born, the CO sublimation radius extends to 100 AU, and the temperature at r < 10 AU is 
higher than 100 K, at which some large organic species start to evaporate. 

We investigated the radial distribution of molecules at 9.3 x 10^ yr after the birth of 
a protostar, when the tempcratiuc is higher than 20 K and 100 K at < 3.9 x 10^ AU and 
< 100 AU, respectively. The time taken by the infalling shells to warm up from 10 K to 100 
K is important for the production of large molecules. We found that this warm-up phase 
occurs for only several 10^ yr in our dynamical model. Once the shells enter the region 
where T > 100 K, they fall into the central protostar in the very short time of ~ 10^ yr. 
Large organic molecules such as CH3OH, HCOOCH3 and CH3OCH3 are formed on grain 
surfaces at temperatures of 20— 40 K and subsequently released into the gas phase by thermal 
evaporation. In our model, only the abundance of HCOOH seems to be infiuenced by gas- 
phase chemistry. As a consequence, the radius at which the gas-phase abundance of an 
organic molecule typically increases strongly in the envelope corresponds to its sublimation 
radius. Our model also indicates that carbon-chain species can be formed by a combination 
of gas- phase reactions and grain-surface reactions after the sublimation of CH4. 

We compared molecular abundances in an isolated core model, which is exposed to the 
ambient UV field, and an embedded core model. We found that the photo-reactions in ice 
mantles are important in determining the ice abundances and thus the molecular abundances 
in hot corinos. The molecule CO2, in particular, is enhanced in the isolated model, while 
methanol and formaldehyde are more abundant in the embedded model. Our model shows 
acceptable agreement with gas-phase observations of the hot corino IRAS 16293-2422 and 
ice-phase observations towards the low-mass protostar Elias 29, considering uncertainties 
and variations in observed abundances. A more detailed comparison, including a radiative 
transfer treatment to directly compare with the observed spectra, is desirable. 
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Table 1: Elemental abundances with respect to H. 



Element 


Abundance 


He 


9.75(-l)^ 


N 


2.47(-5) 





1.80(-4) 


c+ 


7.86(-5) 


s+ 


9.14(-8) 


Si+ 


2.74(-9) 


Fe+ 


2.74(-9) 


Na+ 


2.25(-9) 


Mg+ 


1.09(-8) 


P+ 


2.16(-10) 


C1+ 


1.00(-9) 



"a(-6) means a x 10 



Table 2: Percentage abundances of ice-mantle species compared with H2O ice in model cores 
at tfinai- For each core model, ice abundance at radius of 1000 AU and 8000 AU are listed. 
Ratios of column densities integrated towards the core center are also hsted as "column". 



Species Embedded core Isolated core 

1000 AU^ 8000 AU column 1000 AU 8000 AU column 



CO 


7 X 10-^^ 


10 


0.4 


8 X 10-*^ 


23 


0.04 


C02 


10 


4 


4.2 


125 


21 


45 


H2C0 


6 


3 


1.5 


2 X 10-^ 


0.08 


0.03 


CH30H 


3 


2 


2.7 


0.1 


0.04 


0.2 


CH4 


3 X 10-9 


23 


1 


2 X 10-9 


5 


0.4 



"radii 
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Table 3: Gas-phase molecular abundances in IRAS 16293-2422 and model results. 



Species 


IRAS 16293-2422 


model 






embedded 


isolated 


H2CO 


1.0(-7)^, l.l(-7)^ 


2.8(-6) 


1.3(-11) 


CH3OH 


1.0(-7)'i, 9.4(-8)^ 


3.0(-6) 


8.5(-8) 


HCOOCH3 


2.5-5.5(-7)^ , 2.6-4.3(-9)^ , > 1.2(-8 


)s 1.8(-9) 


3.2(-ll) 


HCOOH 


6.2(-8)% 2.5(-9)s 


1.7(-8) 


2.1(-8) 


CH3OCH3 


2.4(-7)^ , 7.6(-8)" 


3.5(-10) 


5.8(-12) 


CH3CN 


1.0(-8)^ , 7.5(-9)^ 


3.0(-8) 


1.3(-8) 



"Gas-phase abundances at r = 30.6 AU are listed, because the abundances are mostly constant at r < 100 
AU. 



Maret et al. I (|2004 



Chandler et al. I (|2005f) 



' iMaret et all |2005l) 
( 2003h 



Cauzax et al. 



iKuan et al. I (|2004h 

feemiiian fc Hollis I |2006l ) 
ISchoier et al~l (l2002ah 
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Fig. 1. — Evolution of a star- forming core 
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Fig. 2. — Radial distributions of density, temperature and infall velocity in the prestellar 
core (left panels) and protostellar envelope (right panels). The times shown are relative to 
the birth of the second hydrostatic core, the protostar. 
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Fig. 3. — Temporal variation of the central temperature of the prestellar core plotted as a 
function of the central density. 
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Fig. 4. — Sublimation radius of CO (r20K) and large organic species (riooK) as a function of 
total luminosity of the core. The solid and dashed lines depict the radius inside of which the 
temperature is higher than 20 K and 100 K, respectively. 
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Fig. 5. — Temporal variation of density and temperature in infalling shells that reach r = 
2.5, 15, 125 and 500 AU at the final stage of our model (tcore = 9.3 x 10^ yr). The horizontal 
axis is the logarithm of tgnai — ^core- 



- 31 - 



— r-i — 1 1 |-i — 1 1 |-i — 1 1 r 



H2O ice 



H2O gas 




lE+5 lE+4 lE+3 lE+2 10 
tfinal - tcore [yr] 



Fig. 6. — Temporal variation of molecular abundances in the shell that reaches r = 2.5 AU 
at tfinal- The horizontal axis is the logarithm of tgnai — ^core- Black lines represent ice mantle 
species, while gray lines represent gas-phase species. 
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Fig. 7. — Radial distribution of (a) physical parameters and (b-d) molecules in the protostel- 
lar core at tgnai- Black lines represent ice- mantle species, while gray lines represent gas-phase 
species. 
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Fig. 8. — Radial distribution of simple gaseous molecules at tfinai- 
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Fig. 9. — Radial distribution of molecules in the protostellar core, as in Figure 7 (b-d). The 
core is assumed to be directly irradiated by the interstellar UV radiation. 
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Fig. 10. — (a-b) Temporal variation of carbon-chain abundances in the shell that reaches 
r = 2.5 AU at tfinai, ^is in Figure 6. (c-d) Radial distribution of carbon-chain species in the 
protostellar core, as in Figure 7. 



